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Abstract 

We study perturbation of holographic dark energy and find it be stable. 
We study the fate of the universe when interacting holographic dark energy is 
present, and discuss a simple phenomenological classification of the interacting 
holographic dark energy models. We also discuss the cosmic coincidence prob- 
lem in the context of holographic dark energy. We find that the coincidence 
problem can not be completely solved by adding an interacting term. Inflation 
may provide a better solution of the coincidence problem. 
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1 Introduction 

The cosmological constant problempQ is a longstanding problem in theoretical 
physics, and has been taken more and more seriously since the discovery of accelerated 
expansion of our universe[2||3][4j. In addition to the problem why the cosmological 
constant is nonvanishing, there is also the "cosmic coincidence problem" [1] . 

Based on the validity of effective quantum field theory, Cohen et al [5J pointed 
out that the quantum zero-point energy of a system should not exceed the mass of a 
black hole of the same size. This observation relates the UV cutoff of a system to its 
IR cutoff. As a cosmological application, Li [6] suggested to choose the future event 
horizon as the IR cut-off, the energy density of vacuum is given by 



dark energy (HDE) and has been studied extensively [TJ. 

As a phenomenological model, the stability of holographic dark energy is an im- 
portant issue and has been investigated by Myung |8J first. Myung [5] assumed that 
holographic dark energy is a usual fluid component. He calculated the square of sound 
speed of holographic dark energy and found it be negative, this leads to an instability 
of perturbation of holographic dark energy. However, holographic dark energy is given 
by the holographic vacuum energy, whose perturbation should be treated globally. A 
calculation of perturbation of holographic dark energy will be presented in sect. 2. It 
is shown that perturbation of holographic dark energy is stable, and we do not need 
to face the negative sound speed square problem. 

The interacting holographic dark energy, assuming that dark energy interacts with 
matter, has become a popular topic recently [9] . With the interacting term, the story 
of holographic dark energy becomes more interesting. There is often an attractor 
solution to the evolution equation, in which the effective equations of state of dark 
energy and matter become identical in the far future. In sect. 3, we will give a simple 
and phenomenological classification of the interacting terms, we will show that we 
can tune the interacting parameter to avoid the phantom-like universe. 

In sect. 4, we shall discuss the coincidence problem of holographic dark energy in 
a more natural way. We end this paper with conclusion and discussion in sect. 5. 
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2 Stability of Holographic Dark Energy 



In this section, we investigate perturbation of holographic dark energy. First, we 
shall calculate perturbation of the future event horizon, consequently we get the den- 
sity perturbation of holographic dark energy. Finally, we couple this density pertur- 
bation to gravity. As an application, we analyze the coupled equation approximately 
in the dark energy dominated era, and show that the perturbation is stable. We also 
solve the perturbation equation numerically outside the horizon. When dark energy 
dominates, the solution agrees with the analytic result. When matter dominates, the 
numerical result also shows that the perturbation is stable. 

We consider the scalar type perturbation of the metric. In the Newtonian gauge, 
the perturbated metric takes the form 

ds 2 = -(1 + 2$(r, t))dt 2 + a(t) 2 (l - 2$(r, t))dx 2 , (2) 

where for simplicity, we have assumed that the perturbation is spherically symmetric, 
$ = $(r, t), where r = |x|. In this metric, light traveling from the horizon towards 
the origin r = still goes straightly. As illustrated in Fig. [TJ the future event horizon 
R h can be written as 

rr h (t) 

R h (0,t) = / a{t){l-${r,t))dr , (3) 
Jo 

where Rh(0,t) denotes the future event horizon at position r = at time t, Th(t) 
denotes the coordinate distance of the future event horizon. At the first order in 
the perturbation theory, r h = r h0 + 5r h , where r h0 (t) = dt' / a(t') , and Srh can be 
written as 

So the variation of the future event horizon Rh at the position r = takes the form 

5R h (0,t) = R h (0,t) - R h0 = a(t) |^°° ^j0^dt' - jP° S(r,t)drJ . (5) 

Note that for the background value, we have Rho = arho- 

Using the definition of holographic dark energy ([1]), and varying R^, we get the 
variation of the energy density of holographic dark energy with respect to the metric 
perturbation 

5p D = -2p D —± . (6) 



Figure 1: This figure illustrates how to calculate the perturbation of the future event 
horizon. We first integrate along the blue line, which is the geodesic for photons, to get the 
coordinate distance of the future event horizon r^, and then integrate along the green line 
from r = to r = to get the physical distance R^. The dashed red and yellow curves 
illustrates the sub-horizon and super-horizon perturbations, respectively. 

Inserting this equation into the 00-component of the perturbated Einstein equa- 
tion, one obtains 

Z!$ _ 3 #$ _ 3^ = ^(5 PD + 5 Pm ) . (7) 

For simplicity, we neglect the matter density perturbation, so 5p m = 0. To solve this 
equation, we expand $ using its eigenf unction. Write 

^ / \ sin(/cr) . , 

S(r,t) ^ , (8) 

k 

where we have dropped the cos(£r)/r terms, which lead to a singularity at r = 0. 
Then satisfies 

Pd r [a 2 

2<$> k {t')sm{kr m {t'))dt' _ pW sin(fcr) ir _ 

One way to deal with this equation is to take derivative with respect to t. This 
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The Lo^arrtlmi of the Scale Factor x=lii(a'| 



The Logaiitlun of the Scale Factor x=liil a | 



Figure 2: We plot the gravitational potential $ as a function of x = In a in the kr^ -C 1 
case. The left figure is plotted starting from the dark energy dominated era £Id\x=o = 0.72, 
and the right figure is plotted starting from the matter dominated era £Id\x=o = 0.01. The 
blue (dashed), red (solid) and yellow (dotted) curves corresponds to c = 0.8, c = 1.0 and 
c = 1.2 respectively. In both figures, we see that the perturbation approaches a constant 
mode. Note that we have chosen the initial condition $(0) = and $'(0) = 1. So although 
in the right figure, the amplitude can grow 0(100) times before approaching the constant 
mode, but as the initial condition should be set to $(0) oc po, the physical amplitude is 
still not too large. 

integral equation becomes a differential equation 

+ 7777 U 3H + 9H 2 - — + — + - X" / r - dr )■ 



3H [\ R h0 a?) M*R h0 J kr 

This equation can be solved at least numerically. As an example, the evolution of the 
krho <C 1 mode is shown in Fig. [2J 

Eq.Q can also be treated directly in the dark energy dominated era po > p m - 
To investigate the stability of the perturbation mode, we focus on behavior of 
When ~~ * 0, the perturbation mode is frozen, and when ^fc/^fc < 0, the 
perturbation mode is decaying. 

Let us first consider the super-horizon mode kr h0 (t) C 1. In this case, Eq. Q) 
can be written as 



3H ^- 3 ^ k -^ t + J ^-J t wftV) . (11) 



The integral in the last term can be estimated using 

r^^r^*^*^-'') • <i2) 

Note that for f f °° ■^^(t'—t), using dt'/a(t') = —dr h0 and expanding R h0 (t') ~ R h0 (t) + 
[if — t)R h0 (t), we have 

f°° dt' . s f°° dt' . x f [K \ 

I Wf ~ t] = ™ Rm + 1 MFf ~ t] Hfo - l ) • (13) 

where p c is the critical density. So we see that for the dark energy dominated era and 
c ~ 1, the expansion works well. Up to the leading order, J t °° ^tj(£' — t) — r h0 R h0 . 
So Eq. (Q can be written as 

(2c A /^-l)$ fe = tf(l-^)$ fc . (14) 

V Pc Pc 

During the dark energy domination, 1 — £2. approaches zero quickly, so the super- 
horizon perturbation approaches a constant. There is no instability for the perturba- 
tion. 

Note that when c ~ 1/2, there exists a parameter region where 2c\J Pd/Pc — 1 
approaches zero before po —> pc- In this case, the next to leading order correction to 
(13) should be considered. However, the experimental data indicates that c does not 
lie in this regime. 

For the sub- horizon mode kr^it) ^> 1, similar analysis can be performed. One 
can divide the nonlocal integration in into two parts, namely krho{t') ^> 1 and 
krhoit') <^ 1, and use | sin[fcr/ l0 (t')]| < 1 and | sin[/cr/ l0 (t')] I < k r ho(t') respectively. It 
can be shown that the dominant contribution comes from 

3H$ k ~-(^+3H 2 )<f> k , (15) 

and the other terms are suppressed by a factor of l/[krho{t)]. So the sub-horizon 
mode is a decaying mode. Again, no instability appears. 

Before proceeding to the next section, we would like to discuss two physical issues 
of the holographic dark energy. 

First, from the calculation above, we see clearly that perturbation of holographic 
dark energy is nonlocal. This is completely different from a usual fluid component. 
For a usual fluid component, the perturbation equation has a non- vanishing local 
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limit. In this limit, the perturbation equation follows from local conservation of the 
energy-momentum tensor d^T^ v = 0, and the sound speed c s for the perturbation is 
given by c 2 s = dp/ dp. When c 2 s = dp/ dp < 0, the perturbation of the fluid is unstable. 
But for holographic dark energy, the perturbation of the energy density comes from 
the perturbation of the metric, and does not suffer such instability. 

Second, although we discussed the evolution of perturbation for holographic dark 
energy, we did not discuss the initial condition for it. As it is difficult to write the 
holographic dark energy component into the Lagrangian, the quantum initial condi- 
tion for holographic dark energy is not available. Another source for perturbation 
of holographic dark energy is perturbation of the matter component. Perturbation 
of the matter component couples to the metric perturbation, thus providing a initial 
condition for the holographic dark energy perturbation. 

3 The Fate of Interacting Holographic Dark En- 
ergy 

In this section we will make a simple and phenomenological classification of inter- 
acting holographic dark energy. We also study the fate of the universe with interacting 
holographic dark energy: in what case it will be phantom-like, in what case phantom 
will be avoided, and whether the big rip will happen or not. 

For simplicity, use w to denote the effective index of the equation of state of dark 
energy (which is sometimes written as tuff in the literature), and use w m to denote 
the effective index of the equation of state of matter (which is sometimes written as 
in the literature) in the following discussion. Please note that we neglect the 
curvature of the universe, the following calculation was done with the assumption 
that the universe is flat. 

3.1 Dark Energy Decay to Matter 

First, we consider the case when holographic dark energy decays to matter 

p' D + 3(1 + w D )p D = 3bp D , (16) 
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p' m + 3p m = -3bp D , (17) 

where the prime denotes derivative with respect to In a. It is worthy to note that the 
total energy is conserved in the interacting holographic dark energy model, although 
dark energy and matter are not conserved separately. For the lack of the first principle 
of holographic dark energy, we take the above equation phenomenologically. We 
simply follow other works done about the interacting holographic dark energy [9]. 
The decay rate is proportional to the energy density of dark energy, so naturally we 
have b < 0, meaning that dark energy decays to matter. Moreover, b > will lead to 
unphysical consequence in physics, such as p m will become negative and Qd will be 
larger than 1 in the future. So we assume b < in this subsection. 
Comparing with the effective equation of state 



p' D + 3(1 +w)p D = , p' m + 3(1 + w m )p m = 

we find the indices of the effective equation of state 

w = w D -b , w m = bQ D /Q m . (19) 

If the index of the effective equation of state of dark energy satisfies w < — 1, dark 
energy is phantom-like. 

Taking derivative of Eq.([T]) with respect to In a, we have 

p> D =2p D (^°-lj , (20) 

from Eqs.(fTO])(j20jh we get 

w D = ----^^ + b . 21 

6 6 C 

Using the definition of and taking derivative of Qd with respect to In a, we have 

H' 

fl' D = -2Q D + — 5- - 2Q D — . (22) 
c H 



From Eq.fT22l. we get 



H 2tt D 
From the Friedmann equation 



h' _ n' D . 1 + v^_ (23) 



H = -4vrG(p +p) = -AnG(p + p D w D + p r w r ) , (24) 
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we get 

the last term of the RHS of Eq. (1241) can be neglected. 

Substituting Eq. (j25p into Eq. fl22l . we obtain the differential equation for 

^■D ( i n \/t , 2 V^D 



(i-n D )(i + ^—^) + sm D . (26) 

We solve Eqs.(25)(26) numerically and the evolution of the universe has been 
shown in Fig. [3l Now we will discuss these equations analytically. 

Considering Eq.( l26l) . we find that the LHS of Eq. (l26l) will vanish only when the 
scale factor goes to infinity. To see this, we define 

f(y) = 2 -^ = (l-y 2 )(l + 2 -^) + 3by\ (27) 
y c 

where y = \/Vt D . The equation f(y) = has three roots. Since /(0) = 1 > and 
/(l) = 36 < 0, there is one root in the region [0,1] at least. We only consider the 
region [0, 1] since it is the physical region, for Friedmann equation in the flat universe, 
the energy density of holographic dark energy and that of matter should be positive. 
So Qd should never go beyond the region [0,1]. We assume that y\ is the first root in 
the physical region [0, 1]. Thus, we find the integral equation 

rVl dv 

^3 7 = lna-lna(e) , (28) 

y\u=\\y-Vi) 

where e is a cut-off at the early universe. We find that the LHS of Eq.( l28|) diverges, 
which means that the scale factor approaches infinity as y — > y\. Combining the 
above equation with the fact that Q' D > when — > 0, we conclude that Q' D will 
remain positive, and approaches zero when the scale factor approaches infinity. 

When the index of the effective equation of state satisfies w = wo — b > —1, i.e. 
< 1, dark energy does not behave like phantom. As we have discussed, is an 
increasing function with the scale factor. So once the no phantom condition < 1 
is satisfied at a — >• oo, it will be satisfied along the whole history of the universe. By 
setting Q' D = in Eq. (l26|) . we get b < 1 — c~ 2 . It is the necessary condition to avoid 
the phantom phase. 

We can also prove that b < 1 — c~ 2 is the sufficient condition of no phantom. To 
see this, we first investigate the limiting case 6 = 1 — c~ 2 , then use the monotonicity 
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Figure 3: This figure illustrates evolution of the universe with assumption that holographic 
dark energy can decay to matter. The red (solid) curve corresponds to the phantom-like 
universe, and blue (dashed) curve corresponds to the other case. We choose initial condition 
^do = 0.73, ho = 75 in our numerical calculation. 

of Qd in 6 to prove for the general case. Substituting 6=1 — c~ 2 into f(y) = 0, we 
get three roots, y\ = c, y 2 ,3 = ~ 3± ^ 9 ~ 8c -. Note that ?/2,3 < 0, so y\ is the only root 
making sense physically. Note that solution yi = c corresponds to w = — 1, which is 
the boundary of the no phantom condition. 

To see the monotonicity, suppose is the root of f(y) — 0, 

(l-^ 2 )(l + ^) + 36y 2 = 0. (29) 

Taking derivative of the above equation with respect to 6, we get, 

<hh = zM. (30) 

db -2^ + 1-^ + 6^' 

There is only one inflexion, -# = 0, yi = 0, so the function varies monotonously when 
yi > 0. Substitute 6=1 — c~ 2 into Eq. (!30l) . one finds -jr > 0, around yi — > c. That 
means y = c is largest root of Eq.( l29l . So we conclude that for general 6 satisfying 
6 < 1 — c~ 2 , we have < 1 ; and phantom will be avoided. Thus, we conclude that 
6 < 1 — c~ 2 is the sufficient and necessary condition of no phantom. 

In dark energy models, phantom usually causes the big rip. In the remainder 
of this subsection, we will verify that this statement is also true for the interacting 
holographic dark energy. First, we will prove that if the no phantom condition 6 < 
1 — c -2 is satisfied, there will be no big rip. 

From equation 

H=-4nG[(l + w)p D + (l + w m )p m ], (31) 
10 



where w m = ^ > -|(1 + ^^-) > -1, we get H < 0, the Hubble parameter will 
become smaller and smaller and big rip will never happen. 

On the other hand, if dark energy is phantom-like, we can show that the big rip 
will definitely happen. Consider the asymptotic behavior of the evolution equation. 
When a — > oo, to have phantom, we have 1 + w = 1 + w m = —a, where a is a 
positive constant. Eq. (!3T!) can be rewritten as, H = — 4nG[~ ot(po + Pm)] — ^-H 2 in 
the a — > oo limit, and H = HQ _\ at / 2 i where H is a integral constant. So the big rip 
happens in a finite time t — f^ 11 - 

3.2 Hybrid Interaction 

For the hybrid interaction, the interacting term is proportional to the critical 
energy density. The evolution equations are 

p' D + 3(1 + w D )p D = 3bp c , (32) 

Pm + 3 Pm = -36p c . (33) 

Dark energy will become dominant in this case, and po —> p c - For the same reason 
as in Subsection 3.1, we consider the case b < 0. 

With a same procedure as in the previous subsection, we obtain differential equa- 
tions for and H 



-Ji = (l-n J ,)(l + ^-^)+36, (34) 
\lu c 

£ = £ + («) 
Solving above two equations numerically, We get the evolution of the universe, which 
has been shown in Fig. HI The indices of the effective equations of state of holographic 
dark energy and matter are 

1 2y/H^ b 

w = ~a > w m = — . (36) 

Using Eq. (134D and the condition w = — | — 2y ^° > —1, we get the sufficient condition 
of no phantom is b < c 2 — 1. This condition is also the necessary condition, as can be 
shown in the same way as in the previous subsection. We rewrite Eq. (1341) in terms of 
y = y/th,, 

/(y)^ = (l-y 2 )(l + ^) + 3&. (37) 

y c 
n 




Figure 4: This figure illustrates evolution of the universe with assumption that the inter- 
acting term is proportional to the critical energy density. The red (solid) curve corresponds 
to the phantom-like universe, and blue (dashed) curve corresponds to the other case. We 
find that evolution is very similar with the case dark energy decay to matter. 

Assume that yi is the root of f(y) = 0, 

(1-^ 2 )(1 + — ) + 36 = 0. (38) 
c 



2 1 „. _ „ „, _ -3c±Vl6~15c 2 
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There are three roots in the limit case b = c — l,yi = c, 1/2,3 

In addition, Eq. (1281) implies that there is a constraint on b. From Eq. (!28|) . we 
know that holographic dark energy evolves monotonously along the whole history of 
the universe. Take into account the fact that Qd is close to zero at the beginning of 
evolution of the universe, Q' D must be positive all the time. To meet this requirement, 
we need b > —1/3. Otherwise, flu will be a decreasing function in time. Combining 
this with Eq.(EED, we obtain \ < 

Taking derivative of Eq. (l38l with respect to b, we get 

dyi -3 



db ~ ~^-2y l + l 



(39) 



We see that the denominator is negative since - < ^r+yi. So yi increases monotonously 

with respect to b in the physical region. If we tune b such that b < c 2 — 1, we have 

1 _ 2 
3 3' 



w = — \ — < — 1 and phantom phase will be avoided. 



By a similar analysis as in the previous subsection, we find that b < c 2 — 1 is also 
the sufficient and necessary condition of no big rip. 
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Figure 5: This figure illustrates evolution of the universe with assumption that matter 
can decay to dark energy. The red (solid) curve corresponds to the phantom-like universe, 
and blue (dashed) curve corresponds to the other case. We find that dark energy will be 
dominant and = 1 eventually. 

3.3 Matter Decay to Dark Energy 

Finally, we consider the case when matter decays to dark energy, 

p' D + 3(1 + w D )p D = 3bp m , (40) 

p' m + 3 Pm = -36p m • (41) 
In this case, b > 0, and the differential equation for takes the form 

^ = (i - n D )(i + ^) + 36(i - a D ) ■ (42) 

The differential equation for H can be written as 

H 2 c 2 V D! 2 y ! 

We obtain the indices of the effective equation of state 

w = -\- 2 -^, w m = b. (44) 

Taking Q' D in Eq. (f42l . there is only one root in the range [0,1], dark energy will be 
dominant, and Qd approaches 1 eventually. The sufficient and necessary condition of 
no phantom is c > 1. It is also the sufficient and necessary condition of no big rip. 
The analysis is just a little different from the previous ones, because w m = b all the 
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time and there is no attractor solution to the evolution equation of state. Consider 
equation, 

H=-AirG[(l + w)p D + (l + b)p m ], (45) 

If no phantom, i.e. 1 + w > 0, H < all the time, and the big rip will never happen. 
If dark energy is phantom-like, c < 1, Note fljj — > 1 eventually, we can neglect the 
second term in Eq. fj4"5"l) . H = u _^-^-i)t when t = c ^i°_ 1 , the big rip will occur. The 
evolution of the universe in this case has been shown in Fig. [5j 

4 The Coincidence Problem 

The usual solution to the coincidence problem is to calculate the ratio of dura- 
tion of coincidence state and lifetime of the universe, assuming that dark energy is 
phantom-like and the universe will end with the big rip [10] . The coincidence problem 
is solvable if this ratio is not too small. Another solution to the coincidence problem 
is the interacting dark energy [9]. We shall show that the interacting dark energy does 
not solve the coincidence problem in this section. We consider the solution proposed 
by Li in the original paper of holographic dark energy [6] a more natural solution. 

We start with evolution equations 

p' D + 3{l + w)p D = , (46) 

p' m + 3(1 + w m )p m = , (47) 

where the prime denotes derivative with respect to In a, and w, w m are the effective 
indices of the equations of state of holographic dark energy and matter, respectively. 
We rewrite Eq.( )46l) in terms of an integral form 

pin a 

\npi)/pDo= / — 3(1 + w)d\na (48) 
Jo 

Recall the median law of integral, we write the median value of w as a constant w 
varying in the interval (— |, — | — ^). So the above integration can be written as 
hip D /p D Q = — 3(1 + w) lna, we get p D = p D0 a~ 3 ( 1+w \ Similarly we write the median 
value of constant w m , and we get p m = p m oa ( + w ™>^ where poo and p m o are 

energy densities at the time where we set the scale factor ao = 1. 
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From the Friedmann equation 

H 2 = —( Pm + PD ) , (49) 

and the definition of the Hubble parameter H = ^jp, the age of the universe can be 
written in an integral form as follows 

f a ° RttC 

t= / [^( Pm + PD)]- 1/2 d\na. (50) 



o 



3 

Rewrite it in terms of r = t -, we obtain 



i+a i fro 1+t j 

t = # _1 (1 + r ) 1/2 rj* s=s * ;n / (1 + r)- 1/2 r"2(»-»™) _1 dr , (51) 

-3(w - w m ) Jo 

The integration in Eq. (l5TT) satisfies the follow inequality 

(l + r ) 7 r 2 c*-«m) < / (l + r)~ 1/2 r 2 (»-«m) dr , (52) 



o 



and 



So we get 







l + r)~ ' r 2 (™-™™) dr < / r 2 ( fi -»m) . (53) 



HfTo-^l + w^- 1 < t < \h \1 + r ) 1/2 (l + w^- 1 . (54) 



If we input the value of ro at the present time j3] [I] , we find that t is about Hq ( 1 + 
w m )~ l . From the previous section, we know that w m is proportional to b, the inter- 
acting parameter. We see the coefficient of Hq 1 is b dependent. For b = 0, t is about 
-f^ 1 . For b 7^ 0, from Eq. (fl9|) . u? m = for dark energy decaying to matter case, 
and from Eq. (l3~6l) . w m = -J*- for hybrid interaction case, where b is always negative. 
From Eq. ()44l) . we have w m = 6 (b > 0) for the case in which matter decaying to dark 
energy. 

The key issue of the coincidence problem is why the ratio of holographic dark 
energy to matter is order one nowadays, in other words, why the ratio is order one 
when the age of universe is about 10 10 years. The order of magnitude of the age of 
the universe is determined by Hq . The relationship between the age of the universe 
and H^ 1 depends on the initial condition of the universe. The interacting term can 
change the coefficient in the front of H^ 1 , but it can not provide any information 
about the value of Hq 1 . Thus, the interacting holographic dark energy can not solve 
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the coincidence problem completely in this sense. The above approach just solves the 
coincidence problem partially. 

To solve the coincidence problem completely, the initial density of holographic 
dark energy and the influence of inflation should be taken into account. As proposed 
in the original paper of holographic dark energy [6J, the initial energy density of 
holographic dark energy has been inflated away by a factor exp(— 2N) in the inflation 
epoch, where N is the e-folding number of inflation. So the ratio between p^ and p r , 
the radiation density, should be about 10~ 52 at the onset of the radiation dominated 
epoch, if we suppose that the inflation energy scale be 10 14 GeV and inflaton energy 
completely decays into radiation at the end of inflation. This will lead to the order 1 
ratio of holographic dark energy to matter in our epoch. Thus inflation not only solves 
the traditional naturalness problems and helps to generate primordial perturbations, 
but also solves the cosmic coincidence problem. 

5 Conclusion and Discussion 

In this paper, we study perturbation of holographic dark energy. Since holographic 
dark energy is just the holographic vacuum energy, its perturbation is global. We 
calculate perturbation of the holographic dark energy, and find it stable. 

Many numerical and analytic works have been done on the interacting holographic 
dark energy. We made a simple and phenomenological classification of interacting 
holographic dark energy in this paper, and derived the sufficient and necessary con- 
dition of no phantom (the big rip). Needless to say, this classification has not been 
done previously. It is worth to note that we write the interacting term just by hand 
for lack of knowledge of the first principle of holographic dark energy. We hope we 
will return to this issue in future projects. 

We also discussed the coincidence problem. It is shown that the interacting holo- 
graphic dark energy approach only solves the coincidence problem partially. The 
original solution to the coincidence problem proposed in [6] stands a better resolu- 
tion. 
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